ELASTIC-PLASTIC  RESPONSE  OF  STRUCTURES 
TO  BLAST  AND  IMPULSE  LOADS 


Joshua  E.  Greenspon,  Dr.  Eng. 


STATEMENT  NO.  1 

Distribution  of  This  Document  Is  Unlimited 


D  P 

-inr 

4  1967 


Ballistic  Research  Laboratories 
Aberdeen  Proving  Ground 
Contract  No.  DA-18-001-AMC-1019  (X) 
Tech.  Rep.  No.  7 
March,  1967 


Reproduction  in  whole  or  in  part  is  permitted 
for  any  purpose  of  the  United  States  Government. 


J  G  Engineering  Research  Associates 


3B31  MENLO  DRIVE 


BALTIMORE  15.  MARYLAND 


This  report  discusses  the  general  types  of  failures  of  typical  structures 
that  are  used  in  aircraft  and  missiles.  The  theories  of  elastic  and  plas¬ 
tic  deformation  of  these  structures  are  presented  and  comparison  with 
experiments  of  simple  structures  is  given.  It  is' found  that  the  side  on 
impulse  can  be  used  together  with  rigid-plastic  theory  of  buckled  and 
collapsed  cylinders  to  predict  plastic  deformation  for  a  wide  range  of 
pulse  durations.  For  beam  type  structures  the  impulse  and  rigid  plastic 
approximations  seem  to  only  hold  for  pulses  of  very  short  duration. 
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I.  Introduction 


There  are  two  central  aspects  involved  in  calculating  the  response  of 
structures  to  blast  or  other  impulsive  loading  which  deforms  the  struc¬ 
ture,  The  first  aspect  is  the  description  of  the  loading  and  the  sec¬ 
ond  is  the  method  of  obtaining  the  elastic  and  plastic  response  once 
the  loading  is  known.  The  classical  method  of  obtaining  the  response 
of  any  structure  is  to  determine  the  differential  equations  governing 
the  behavior  of  the  structure  and  solve  these  by  some  technique  such 
as  a  finite  difference  procedure,  the  Galerkin  method,  or  obtain  exact 
solutions  if  possible.  Such  classical  solutions  are  based  upon  the 
premise  that  the  loading  on  the  structure  is  known  completely,  i.e. 
that  its  amplitude  and  its  spatial  and  temporal  characteristics  are 
known.  There  is  much  work  to  be  done  in  determining  both  the  loading 
and  the  response  to  a  completely  defined  load;  especially  to  the  point 
where  the  analysis  can  be  used  with  complete  confidence  in  practical 
situations.  It  is  the  purpose  of  this  report  to  review  some  of  the 
recent  procedures  that  have  been  suggested  in  the  past  several  years 
for  computing  the  loading  and  response  of  beam  type,  plate  type,  and 
shell  type  structures  and  to  point  out  some  of  the  advantages  and  short 
comings  by  comparing  the  results  with  some  available  experimental  re¬ 
sults. 

II.  Damage  in  structures 

A.  General  considerations 

Often  in  the  field  of  blast  damage  the  terms  critical  impulse  or  cri 
tical  pressure  are  used  to  denote  values  for  which  the  damage  reache 
a  certain  approximate  value.  The  concept  of  critical  values  is  ex¬ 
tremely  valuable  in  experimental  studies  *such  as  those  conducted  on 
full  scale  structures  or  those  conducted  by  Schuman’'"' 3*  on  model 
shells.  However  one  can  and  should  expect  more  from  a  theory  than 
just  to  predict  critical  values.  In  fact  due  to  the  inherent  non- 
linearities  in  structures  the  answers  that  one  obtains  in  the  large 
deflection  region  is  a  dependency  between  load  (or  some  function  of 
load)  and  deformation. 

B.  Types  of  failures 

1.  Beam  type  structures  (wings,  control  surfaces,  slender  missiles, 
slender  fuselages) 

In  structures  which  can  be  replaced  by  a  line  with  appropriately 
distributed  mass  and  inertia  properties  the  elastic  and  plastic 
behavior  is  rather  straightforward.  The  structure  behaves  elas- 


*Superscripts  refer  to  references  listed  at  the  end  of  the  report 


tically  and  responds  in  its  modes  until  the  yield  point  is  reached  at 
some  section  of  the  beam.  At  this  Section  a  hinge  forms  and  the  struc¬ 
ture  deforms  plastically  around  this  hinge  as  shown  in  Figure  1. 


la»  Elastic  Response  of 
a  Cantilever 


lb.  Elastic  Response  of  a 
Free-Free  Beam 


lc.  Plastic  Response  of  Id.  Plastic  Response  of  a 

a  Cantilever  Free-Free  Beam 

Fig.  1  Typical  Response  of  a  Beam  Type  Structure 

Since  these  problems  are  one  dimensional  and  usually  involve  only  one 
plastic  hinge  they  can  be  solved  by  starting  with  the  differential 
equations  of  an  elastic  beam  in  bending.  An  investigation  was  under¬ 
taken  several  years  ago^  to  study  such  problems  and  a  general  computer 
program  was  developed  to  solve  this  general  type  of  problem.  This  area 
has  been  under  investigation  for  a  number  of  years  and  there  have  been 
many  papers  written  on  the  subject.  From  a  practical  standpoint  this 
type  of  problem  is  fairly  well  in  hand  at  the  present  time  if  the  load¬ 
ing  on  the  structure  is  adequately  defined. 

2.  Plate  type  structures  (Parts  of  fuselages,  plate  type  control  surfaces, 
plating  between  longitudinal  and  transverse  stiffeners) 

The  flat  elastic  surface  has  one  more  spatial  dimension  than  the  beam 
but  its  elastic  and  plastic  deformation  characteristics  are  many  times 
more  complicated.  In  the  elastic  region  the  governing  differential  equa¬ 
tion  usually  has  to  be  solved  by  a  two  dimensional  finite  difference 
network  or  by  some  approximate  modal  method  since  the  exact  solution 
can  be  obtained  only  for  very  few  boundary  conditions.  The  classical 
plasticity  problem,  especially  for  the  rectangular  plate,  offers  an  ex¬ 
tremely  difficult  problem  even  for  very  small  plastic  deformations. 

For  that  reason  approximate  approaches  were  developed  several  years  ago  ' 
for  small  and  large  elastic  and  plastic  deformations  of  plates. 


Cylindrical  Shells 

The  elastic  analysis  of  cylindrical  shells  has  been  the  subject  of  many 
studies  in  recent  years.  However  the  only  plastic  problems  associated 
with  cylindrical  shells  that  have  been  approached  from  the  classical 
viewpoint  of  flow  plasticity  is  the  axially  symmetric  problem.  Other 
problems  have  been  considered  from  the  standpoint  of  deformation  theory. 
The  real  problems  in  blast  deformation  of  cylindrical  conical  and  other 
types  of  shells  or  composite  bodies  are  not  axisymmetric  and  involve  sit¬ 
uations  where  the  pressure  time  relation  on  the  shell  is  not  easily  ob¬ 
tainable.  This  is  the  reason  whv  a  work-energy  approach  has  been  follow¬ 
ed  in  several  recent  reports.^'  In  that  work  the  asymmetry  of  the 
load  and  deformation  was  considered  as  well  as  strain  hardening,  large 
nonlinear  deformations,  and  thermal  aspects.  It  was  shown  experimentally 
by  Schuman^'^'^  that  cylindrical  shells  under  a  blast  which  engulfs  the 
structure  either  fails  in  a  well  defined  buckling  pattern  or  a  collapse 
pattern.  Moreover,  for  loading  which  is  more  of  a  local  type  such  as 
sprayed  explosive  the  deformation  was  almost  uniform  along  the  length 
but  changed  rather  abruptly  to  zero  at  the  supported  edges.  In  addi¬ 
tion  if  the  explosive  was  sprayed  over  half  the  periphery,  then  only 
this  much  deformed.  A  representative  pattern  is  shown  in  Fig.  2. 


Fig.  2  Typical  Deformation  Pattern  from  Sprayed  Explosive 
Over  Half  of  a  Fixed-Fixed  Shell 

Deformation  of  composite  structures 

Actual  aircraft  and  missiles  are  composite  structures  which  can  be  con¬ 
sidered,  in  an  approximate  sense,  to  be  made  up  of  beam,  plate,  and  shell 
components.  In  such  structures  several  types  of  failure  may  occur  si¬ 
multaneously,  and  for  this  reason  careful  judgement  should  be  exercised 
in  evaluating  the  blast  resistant  characteristics  of  such  built-up  struc¬ 
tures.  For  example,  an  aircraft  is  composed  schematically  of  a  beam 
type  fuselage  which  is  built-up  of  a  shell  with  ring  and  longitudinal 
stiffeners.  The  wings  can  be  considered  as  beams  with  variable  mass  and 
stiffness.  The  fuselage  can  be  considered  as  a  shell  or  a  beam  when  ex- 


posed  to  an  engulfing  blast  which  loads  the  entire  structure  but  it  can 
also  be  doiisidered  as  a  series  of  stiffened  and  unstiffened  panels.  By 
the  same  token,  the  wings,  elevator  and  stabilizer  act  either  as  beams 
or  a  series  of  panels. 

In  general  we  can  divide  the  deformation  into  primary,  secondary  and  ter¬ 
tiary  components  as  was  done  for  elastic  deflection  of  ships  a  number  of 
years  ago.^  The  primary  deformation  consists  of  deflection  of  the  entire 
structure  (wing  or  fuselage)  as  a  beam.  The  secondary  deformation  invol¬ 
ves  deflection  of  entire  stiffened  plate  or  shell  sections  and  the  ter¬ 
tiary  deformation  involves  deflection  of  plating  between  stiffeners  or 
shell  deformation  between  rings. 

It  is  important  to  look  critically  at  each  type  of  failure  to  see  if  it 
would  kill  the  utility  of  the  structure  in  performing  its  mission.  Re¬ 
ferring  to  Fig.  3  it  is  seen  that  deformation  can  take  place  in  the  follow 
ing  ways: 

1.  Overall  hinging  of  the  wing,  rudder  and  stabilizer 

2.  Deformation  of  the  panels  contained  in  the  fuselage,  wing  and  tail 
assembly 

3.  Buckling  or  failure  of  a  section  of  the  fuselage  and  consequent  hing¬ 
ing  of  the  fuselage  around  the  weakened  buckled  section 

In  a  missile  type  structure  shown  in  the  same  figure  the  deformation 
patterns  are  similar. 

Lumped  parameter  vs.  classical  theories 

In  classical  theory  the  differential  equations  of  the  structure  are  set 
up  for  the  elastic  and  plastic  ranges  and  the  appropriate  pressure  dis¬ 
tribution  in  space  and  time  is  ulped.  The  resulting  answers  are  deflec¬ 
tions,  stresses,  and  strains  as  a  function  of  time  at  many  locations  on 
the  structure.  This  type  of  theory  for  predicting  elastic  and  plastic 
behavior  is  practical  for  beam  type  structures  where  only  one  deflection 
and  only  one  space  dimension  is  present.  For  plate  and  shell  structures 
where  two  space  dimensions  are  present  and  more  than  one  deflection  may 
be  necessary  to  describe  the  behavior  of  the  structure,  the  analysis  is 
many  times  more  involved.  For  many  cases  such  as  large  lateral  deforma¬ 
tions  of  plates  and  cylindrical  shells  an  equivalent  energy  procedure  and 
variational,  approach  such  as  used  in  previous  reports^  offers  many  sim¬ 
plifications.  In  the  equivalent  energy  approach  the  estimated  energy  de¬ 
livered  from  the  explosion  is  equated  to  the  internal  work  done  by  the 
structure  in  plastically  deforming.  The  main  unknown  which  has  to  be  in¬ 
serted  into  the  theory  is  the  deflection  pattern  in  the  plastic  region, 
which  is  determined  from  tests.  Once  this  deformation  shape  is  defined, 
it  is  straight  forward  to  compute  the  work  done  by  internal  forces  in 
the  structure  to  deform  plastically  in  this  pattern.  In  many  cases  it  is 
the  maximum  damage  or  plastic  deformation  which  is  of  primary  significance 
The  lumped  parameter  or  energy  type  analysis  is  capable  of  predicting  the 
necessary  informaton  for  many  situations. 
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Fig.  3  Schematic  of  Various  Types  of  Failure 


It  will  be  seen  later  that  for  short  pulses  the  above  impulse  or  energy 
approach  is  accurate,  but  for  some  cases  involving  longer  duration  blasts 
we  have  to  make  recourse  to  a  more  accurate  procedure  which  uses  the  spa¬ 
tial  and  time  characteristics  of  the  load.  The  analysis  which  fits  this 
qualification  is  the  variational  approach.  In  this  approach  a  variational 
principle  associated  with  deformation  tVpe  (as  opposed  to  flow  type)  plas¬ 
ticity  is  used;  the  principle  is  essentially  Hamilton's  principle  inter¬ 
preted  for  the  plastic  region. 


IV.  Fundamental  energy-impulse  relations  and  the  variational  principle 


A.  Energy  -  impulse 


Let  V  represent  the  work  done  by  the  internal  forces  while  the  struc¬ 
ture  is  deforming  and  let  E  represent  the  total  energy  delivered  to 
the  structure  from  the  explosion;  then  the  fundamental  energy  relation 
that  will  be  used  is  as  follows; 

For  cases  where  the  impulse  is  given  or  can  be  obtained  more  easily 
than  the  energy,  an  alternate  relation  is  derived  below. 


Let  I  be  the  impulse  per  unit  mass  applied  to  the  structure.  The  im¬ 
pulse  momentum  relation  for  an  elemental  mass  dm  can  be  written  (ne¬ 
glecting  the  tangential  and  longitudinal  velocity,  assuming  that  the 
lateral  velocity,  ^ r  is  much  greater  than  the  other  two) 


Thus 


r 


=  x 


[2] 


The  Kinetic  energy  imported  to  the  structure  is 

is'*.?)  *  i  XVA  [3) 

where  ^ )  is  the  mass  per  unit  area  of  the  structure  (which  can 

vary)  and  dA  is  an  elemental  area.  The  impulse  can  vary  over  the  sur¬ 
face  so  write  it  as  follows; 

jzrrx.pj  =  -zzS/%?)  [4] 

Thus 

T~  h. ?)  Z' V*  y.)ctA  [5] 

Equating  this  initial  Kinetic  energy  to  the  energy  of  deformation  ab¬ 
sorbed  by  the  structure  the  expression  for  the  impulse  per  unit  mass 
becomes 


2.  V 


[6] 
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The  total  impulse  on  the  structure  will  then  be 


2  V 


The  impulse  per  unit  area  will  be 


2.  \/ 


If  both  the  mass  per  unit  area  and  the  impulse  are  uniform  then 

**•_  \  / 


*=  V  ^r- 

where  A  is  the  total  loaded  area  of  the  structure. 

B.  The  variational  approach 

The  principle  of  extremum  potential  energy  can  be  applied  to  the  plastic 
regime  if  we  are  dealing  with  a  plastic  deformation  theory  (i.e.  not  a 
flow  theory)  as  we  are  in  this  work.  Greenberg^2  states  this  principle 
as  follows: 

S  fuj  = O 


where 


U  -  J ,  C  f  <n^<rL  )  dv  '  -  f  TiU;  c/ 2  ~  V-  W 

v'  ■4- 


The  integral  over  y‘  is  the  potential  energy  or  work  done  by  the  internal 
forces  during  deformation  and  the  integral  over  51  is  the  work  done  by 
the  external  surface  forces  which  give  rise  to  the  internal  forces. 4 

The  dynamic  counterpart  of  the  principle  of  extremum  potential  energy  is 
Hamilton's  Principle  which  can  be  stated  as  follows: 

sf  £ r-uj  ott  =  o  [U 

where  T  is  the  Kinetic  energy  and  V  is  the  potential  energy  minus  the  work 
done  by  the  external  forces . 

We  are  dealing  with  large  lateral  deformations  of  thin  walled  structures 
where  we  know  the  deformation  pattern  from  test  results  but  the  amplitude 
of  lateral  deformation,  w,  is  inknown.  If  we  neglect  the  longitudinal 
and  tangential  inertia  terms  the  Kinetic  energy  of  the  structure  becomes 

r  -  {  £s«c/n^cJA  [12 

ySA  =  mass  per  unit  area  of  the  structure 

c/A  =  element  of  surface  area 

-w*  =  velocity  normal  to  surface  of  structure 

=  f-i JJfaJwhere  denotes  the  time  dependency  of  the  deflec¬ 

tion  and  fSA)  denotes  the  space  dependency 

13 

The  variation  is  taken  just  as  in  the  elastic  problem  as  given  in  Love. 
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VI 


This  problem  is  equivalent  to  the  following  problem  in  the  calculus  of 
variations: 

Find  the  function  y(x)  which  takes  on  a  given  value  for  a  and  b 

and  which  minimizes  the  definite  integral 

x'=  /  pOJ*. 

The  result  is  that  F  must  satisfy  the  Euler  equation 


In  jur  case 


4.(1$)-° 


f=  =  T-U 

ct)  J  X-  -*  t 


-9  /WJ  ) 

Thus  we  have 

1~  -  ■£  AI-a^^A  s  ^  A) 4^9*74)  -f 

In  general  we  can  fit  a  power  series  to  the  potential  energy  (or  work 
done  by  the  internal  forces) 
i.e . 


[13] 

[14] 

[15] 

[16] 


V  =  A  4-  &  <uJ-m  +  £Vt *s-4  2  +  J/  -  -  - 

i?  CAj  *  ) 

where  is  the  pressure  applied  to  the  surface  of  the  struc¬ 

ture 


[17] 

[18] 


The  result  is  the  following  nonlinear  differential  equation  in  time  for 
the  unknown  A^ya 

*J~C  f  a<CA)  J  r  / ?>&,*)  SSX)  J A 

A  vi 

The  initial  conditions  are 

/^~0  fc )  —  'UT,_  j  ~  /^re 


[19] 


[20] 


where  A^c  and  are  the  deflection  and  velocity  at  the  beginning  of 

the  plastic  regime  which  are  determined  from  the  elastic  analysis  when 
the  yield  point  is  reached. 


Instead  of  the  deflection  being  a  function  of  one  particular  shape 
supposing 

✓vu r  -  2* 


then 


T  *  i  J^A*fA)*r-4fA  =  i fAAAfA)[%»-o„(0^fA)]^dA 

V  =  V  C *^0  *w'  >  ) 


[21] 


[22] 


V. 

M 

1 


Unless  we  assume  complete  decoupling  among  the  terms  the  problem  will  get 
completely  out  of  hand;  thus  we  will  assume  that  this  problem  reduces  to 
finding  the  set  of  a°^^s  satisfying 
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fy,X(A)dA  +  J  r  -t )  ^v>  f'A)  c/A  [23] 

with  initial  conditions  on  displacement  and  velocity  determined  from  the 

end  of  the  elastic  region. 

V.  Description  of  the  parameters  associated  with  the  loading 

A.  General  problem 

Since  the  structure  reacts  wifh  the  blast  it  is  usually  very  difficult 
to  determine  the  complete  spatial  and  temporal  characteristics  of  the 
pressure  which  acts  on  the  structure.  When  the  blast  comes  in  contact 
with  the  structure  it  diffracts  (bends)  around  the  structure  in  a  very 
short  time  and  then  produces  a  drag  type  loading  which  lasts  a  much 
longer  period.  In  general  even  the  deformation  of  the  structure  will 
effect  the  pressure  distribution  and  in  many  cases  involving  control 
surfaces  subject  to  Mast  this  is  considered.  In  many  of  the  more 
complicated  case*  it  is  possible  to  use  some  parameter  to  represent 
the  load  other  than  the  pressure  distribution  itself,  such  as  energy 
or  impulse.  In  such  cases  the  structural  theory  must  be  altered  accord¬ 
ingly  in  order  to  compute  characteristics  of  the  system,  such  as  energy 
absorption,  which  are  most  easily  integrated  with  this  load  parameter. 

B.  Pressure  distributions 


There  are  a  limited  number  of  cases  where  pressure  as  a  function  of 
ppace  and  time  has  been  obtained  for  rigid  bodies  of  various  shape. 

The  Nuclear  Effects  Handbook14  gives  a  particularly  good  outline  of 
the  results  to  date.  A  recent  book  by  Kornhauser1 •  also  gives  pressure 
curves  for  some  particular  cases.  For  aircraft  and  missile  structures 
we  can  approximate  the  pressure  distribution  and  time  history  by  re¬ 
presenting  the  shapes  of  the  various  components  of  the  structure  by 
idealized  shapes  such  as  long  slender  boxes  and  cylinders  and  then  use 
the  pressure  distribution  for  these  shapes.  The  wings  and  other  con¬ 
trol  surfaces  of  aircraft  and  missiles  can  be  represented  by  flat  beam 
type  surfaces  (long  slender  boxes)  and  the  fuselages  can  be  represented 
by  cylindrical  shapes.  The  problem  is  one  of  finding  the  loading  while 
the  weapons  are  parked  and  when  they  are  in  the  air.  The  pressure  time 
hij^ory  for  a  blast  on  a  parked  vehicle  is  shown  schematically  in  Fig. 
4. 
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The  various  parameters  associated  with  each  of  these  plots  can  be  obtain' 

ed  from  the  above  mentioned  reference^  and  work  by  Baker  et.  al.  and 
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Baker  and  Schuman.  Extensive  data  on  incident  and  reflected  pressure 
and  time  duration  of  pressures  associated  with  explosions  of  pentolite 
is  contained  in  reports  by  Baker,  et.  al.  and  Goodman. The  use  of 
this  data  will  be  illustrated  later  in  the  report  when  some  specific 
problems  are  considered. 


During  flight  the  problem  has  to  be  handled  as  a  lifting  body  gust  prob¬ 
lem  and  pressure  time  histories  would  look  schematically  as  shown 

in  Fig.  51  '  for  lifting  surfaces.  Although  there  is  some  lift  on  the 
fuselage  it  is  believed  that  the  parked  condition  can  be  used  as  a  first 
approximation . 

C.  Energy  and  impulse  of  an  engulfing  explosion 


1.  Energy 


For  blasts  which  are  at  some  distance  from  the  structure  the  blast 
wave  engulfs  the  entire  target  when  it  reaches  it.  The  characteris¬ 
tics  of  the  blast  are  described  by  a  charge  weight  (or  weight  of  ex¬ 
plosive)  and  a  distance  from  the  target.  The  weights  for  nuclear  ex¬ 
plosions  are  given  in  terms  of  equivalent  weight  of  TNT.^4  For  ex¬ 
plosions  involving  pentolite  instead  of  TNT  the  conversion  is  that 
1.1  pound  of  TNT  is  equivalent  to  1  pound  of  pentolite  and  the  weight 
of  pentolite  in  free  air  is  equivalent  to  1.8  times  its  weights  on 
the  ground.  7  For  underwater  explosions  experiments  have  shown2  that 
the  energy  flux  density  (energy  per  unit  area)  and  the  impulse  per 
unit  area  before  the  shock  wave  hits  the  structure  can  be  written 
functionally  as  follows: 


x  s  cw\  ,  a,,* wi(^y 


[24] 


where  C,  k,  and  &  are  constants  which  are  determined  by  fitting  ex¬ 
perimental  results  to  the  above  relations. 2^  Using  the  work  of  Baker 
and  Schuman^7  a  similar  analysis  can  be  made  on  air  blast  of  pentolite 
explosive  as  will  be  shown  below. 


22 

Following  the  reasoning  of  Keil,  if  it  _s  assumed  that  the  blast 
wave  is  plane  then  the  energy  flux  (energy  per  unit  area)  can  be  writ 
ten 


[25] 


where  J>4  is  the  air  density,  CQ  is  the  sound  velocity  in  the  air  and 
p(t)  is  the  pressure  at  any  given  point  as  a  function  of  time.  If  it 
is  further  assumed  that  the  wave  is  exponentially  decaying,  then 

«  ?>.e~ [26] 

where  pQ  is  the  pressure  amplitude  at  t  =  0  and  (9  is  the  time  con¬ 
stant  of  the  decay.  Using  the  above  two  relations  it  is  found  that 

£  -  [27] 

*  f.C> 
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The  experimental  curves  of  Baker  and  Schuman  can  then  be  lused  to  obtain 
pQ  and  &  as  a  function  of  charge  weight  and  distance.  Substituting  these 
extrapolated  values  of  pQ  and  0  and  using  the  appropriate  value  for  /?, 
and  C0  the  following  final  relation  is  found  for  the  energy  flux  (energy 
per  unit  area)  in  the  blast  field  at  a  distance  R  from  the  center  of  the 
explosion : 


Ef  'X  t  2.  y/O3  ( 


■&  /  >*- 


where  R  must  be  given  in  feet  and  W  in  &  .  It  is  unfortunate  that  we 
have  to  attach  specific  dimensions  to  the  input  and  that  we  have  to  mix 
inches  and  feet  in  the  same  formula.  However,  it  is  most  convenient  to 
present  the  relation  in  this  form  since  it  is  empirical  and  was  obtained 
by  fitting  the  relation  to  data  that  was  given  originally  in  these  spec¬ 
ific  dimensions. 

2.  Impulse 

Values  of  side  on  (i.e.  incident)  impulse  and  reflected  impulse  (i.e.  re¬ 
flected  from  a  flat  rigid  target  of  infinite  extent)  were  studied  exten¬ 
sively  by  Goodman^-8  for  free  air  blast  of  bare  spherical  pentolite. 
Goodman's  formulas  are  as  follows: 


Side  on 


tmi  _  *^3  7  rsi  „,//,. 

L  x  ^  v  j a  J  if* 


[29] 


*  y 


Reflected:  JT__  _  q/p.2.  7^cST  /£o/  n  /osj  ~,t/.sr c  C303 

where  X  is  the  impulse  per  unit  area  (in  psi  milliseconds)  pQ  is  the  am¬ 
bient  sea  level  pressure  in  atmospheres,  W  is  the  charge  weight  in  pounds 
and  X  is  a  dimensionless  number  equal  to 

X-  */,323  IV*  [31) 

in  which  R  is  the  distance  from  the  center  of  the  explosion  to  the  target. 

Impulse  and  energy  content  for  local  explosion  (such  as  sprayed  explosive) 

As  was  seen  above,  for  explosions  that  are  at  some  distance  from  the  target 
the  weight  of  explosive  and  distance  away  determines  the  energy  content  at 
any  point  in  the  field.  For  impulsive  loading  or  any  type  of  explosion 
which  is  an  integral  part  of  the  structure  such  as  sprayed  explosive,  the 
problem  is  not  so  straight  forward.  In  these  cases  it  will  be  assumed  that 
the  impulse  per  unit  area,  applied  to  the  structure  can  be  measured. 

Then  the  impulse  relation  in  Section  IVA  can  be  employed  to  determine  the 
deformation,  i.e.  r—  - 


[32] 
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in  which yM  is  the  mass  per  unit  area,  V  is  the  energy  absorbed  and  frxjf.) 
describes  the  impulse  distribution.  If  the  impulse  and  mass  per  unit 
area  are  both  uniform  then  the  energy  flux  will  be 


Jl  -  y  -  1  (M* 

A  ~  A  ~ 


[33] 


VI.  Work  done  by  internal  forces  and  coupling  with  the  variational  principle 


A.  Beam  type  structures 


1.  General 


The  complete  theory  for  hinging  deformation  of  beam  type  lifting 
surfaces  was  derived  by  the  MIT  group^'^**  and  by  the  present  auth-*. 
or.4'6  Xn  the  event  that  complete  load-time  histories  are  available 
it  is  certainly  preferable  to  use  the  complete  beam  type  theories 
mentioned  above.  However  for  rough  calculations  and  in  cases  where 
overall  damage  criteria  are  being  sought  a  theory  consistent  with 
the  impulse  or  energy  approach  and  the  variational  approach  discussed 
previously  can  be  employed. 

2.  Elastic  region 


Before  plastic  hinging  occurs,  it  is  assumed  that  the  beam  vibrates 
in  its  elastic  modes.  The  kinetic  and  potehtial  energy  of  the  struc¬ 
ture  in  the  elastic  region  can  then  be  written  as  follows: 

KE  t  ~  f 7  j  4j 

PE  V  =  i  /'  '‘f-X./x j  [  ^ 

Assuming  a  deflection  of  the  form 

✓W  -  <2T  )  [351 

we  then  have  (noting  that  cross  product  terms  disappear  due  to  or¬ 
thogonality  of  modes) 


v-  Z 


[36] 


Substituting  into  Hamilton's  principle  (or  into  Lagrange's  equations) 
we  obtain  the  following  differential  equation  for 


] 


where  yUA*.)  is  the  mass  per  unit  length,  jOfiCj't)  is  the  pressure, 

6  (x.)  is  the  local  width  of  the  beam.  The  initial  conditions 
under  the  assumption  of  blast  loading  are 

S' o)  ~  {*)  =  O  [38] 

3.  Plastic  region 

In  the  plastic  region  the  potential  energy  is  replaced  by  the  work 
done  at  the  hinge  as  follows: 
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Thus  we  have  linear  differential  equations  in  both  the  elastic  and  plas 
tic  regions.  The  one  in  the  plastic  region  has  to  be  solved  piecewise 
for  each  segment. 


An  altermate  method  of  dealing  with  the  plastic  region  is  to  represent 
the  M-  O  curve  by  a  power  series  in  &  as  follows: 

A7  —  Cc  -t-C,  &  -f  - *  ' 


[44] 


The  resulting  differential  equation  is  then  the  following  nonlinear  equa¬ 
tion  q 

The  initial  conditions  for  the  plastic  region  are  the  displacement  and 
velocity  at  the  end  of  the  elastic  region.  If  we  use  the  linear  seg¬ 
mented  curve  then  the  initial  conditions  for  each  segment  are  obtained 
from  the  final  conditions  of  the  previous  segment.  Using  the  Allen  in¬ 
terpretation  of  the  "frozen  hinge" ^5  it  is  assumed  that  the  final  plastic 
deformation  of  the  structure  consists  of  a  rotation  around  the  hinge  and 
occurs  at  the  time  when  the  angular  velocity  reaches  zero  for  the  first 
time. 


B.  Plate  type  structures 
1.  Elastic  deformation 


For  unstiffened  or  stiffened  and  sandwich  plates  the  elastic  behavior 
of  the  plate  can  be  computed  by  using  the  theory  developed  in  a  paper 
by  the  author,  published  several  years  ago.-*  The  basic  equations  go¬ 
verning  the  behavior  of  the  plate  is  the  orthotropic  plate  equation 

&  v*  &  *  p**** 


Ot  +  2.H 


a. 


where 


[46] 


In  the  above  mentioned  paper  the  solution  is  given  in  detail  for  a 
rectangular  plate  under  dynamic  loading.  This  elastic  solution  can 
be  written  as  follows: 

*  * 

^  f0  *1 

rt)  =  TV  /  a  -  *)  J-c 

ft.,  =  /  £  4** 

1  *  *  * 


[47] 


[48] 


[49] 


[50] 
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In  the  above  equations  0*.^  P^  and  H  are  the  orthotropic  constants 
which  can  be  determined  from  the  analysis  given  in  an  earlier  paper, 
si*,  is  the  mass  per  unit  area  of  the  plate,  is  the  natural  frequen¬ 
cy  of  the  mode,  %l  and  are  beam  functions  describing  th*  modal 

pattern  in  the  x  and  y  directions  respectively,  fa  and  are  constants 
associated  with  these  mode  patterns,  a  is  the  shorter  side  of  the  plate 
(parallel  to  the  x  axis),  b  is  the.^pnger  side  of  the  plate  (parallel  to 
the  y  axis)  ,  )  describes  the  spatial  distribution  of  the  load  and 

f(t)  describes  the  timewise  distribution. 

2.  Plastic  deformation 

The  theory  of  plastic  deformation  of  isotropic  and  stiffened  plates  is 
given  in  a  recent  report.  The  work  of  deformation,  V  per  unit  volume 
of  an  elastic-plastic  body  can  be  written 

V  =  4-  [51] 

where  <T;  =  V  fqj -ft) V fa-rd  V  T*S) 

<9V  =  -h  Gy.  +  £  i 

The  relation  between  G2  and  <Zi  describes  the  stress  strain  law  of  the 
material.  Assuming  an  incompressible  material  (  &  -  O  )  and  consider¬ 

ing  only  plate  stresses, 

v  L-C  l  1  [52] 

where  =  element  of  volume 


07  =  V<a-<£0>  +<r^+  jy  [53] 

*  =  \Z7777T6~777,  i> ^ 

y 

If  we  consider  that  only  the  lateral  deformation,  'US,  i“s>  of  significance 
and  neglect  the  longitudinal  deflection  u  and  the  transverse  deflection^ 
then  the  strains  ^  and  become28 

,  /  /jw \  *-  _  t)  _  /  /Jur\  **  -  0  ii~r- 

**■  -K ttJ  -  ^-ir  sjy [54] 

^  d  W  2  -  d  A*J~. 

*3-  -  7Z  S%.  ~  2  *  SCSP 

where  £■  is  the  distance  of  any  element  of  the  plate  from  the  neutral 
plane  as  shown  in  the  figure  betow: 


V':' 
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Further  restrict  the  material  to  one  that  obeys  an  elastic-linear  harden¬ 
ing  law  as  shown  below 


Fig.  8  Linear  Hardening  Stress-Strain  Curve 
The  stresses  can  be  written  in  terms  of  strains  as  follows: 


*?■ x  1  '  i  **  J 

•4  fe  *■ 


[55] 


~  "3  77  ‘*r 

in  which  the  stress-strain  law  can  be  written 

Q  *  E  p  -  [56] 

where  s  o  ^ */«*’£•*  ) 

^  0  -  e*/ti )  {ft+'J-'c )  =  /-  g 

Substituting  the  linear  hardening  law  into  the  Expressions  for  the  stress¬ 
es  and  then  into  the  relation  for  V,  we  obtain  the  work  done  by  the  in¬ 
ternal  forces  on  the  plating  of  the  stiffened  plate 


y  *J.J.  J  ,[*£  (•-*)  + 

s ***  FA*  v  v 

where  £  is  the  distance  from  the  neutral  plane  of  the  stiffened  plate 
to  the  midBurface  of  the  fats,  h  is  the  thickness  of  the  plate  and  a,b 
and  the  width  and  length  of  the  plate  respectively.  If  there  is  a  top 
plate  then  there  is  a  similar  expression  V2  that  has  to  be  added  for  this 
part,  i.e. 


[57] 


d  yt/i  t  •  C  S  Sf  c  »*w  ^  V? 


[58] 
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Now  substituting  the  expression  for 

V,  /jffi-*)  '[ 

/■£*/b  fi+4i 


,I?9] 


Substituting  the  expressions  for  the  strains  in  terms  of  deflections, 
we  have  __ ,  V* 

W*/’  -  .  ,  \  *= 


Ji  2'  i 


V-  *2  gj  /TV* k+»  \  A  /  *>  y  a  s 

/T  / 


±£-^*~lY2£±± 


.**  5>t 


-  XT 

and  a  similar  expression  for  V2  if  two  surface  plates  exist 


J  ~  j. 
2  =  * ^ 


In  the  above  equation 


*  =  <  rt&Y&F'-  if  if) 


f  •  (wvmp)>(f.r>f£%  r 


For  the  stiffeners  we  use  an  analysis  similar  to  that  employed  in  ano 
ther  earlier  reference” 
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Fig.  10  Location  of  Stiffeners 
For  a  stiffener  located  at  y  =  yr 

V£.  -  J ' f*  -X  )  £>.  +  E  A  €.  s  ~J  c/Ar  ^ g/fi  f  d*. 


[62] 


Where  A^  indicates  integration  over  the  cross  sectional  area  of  the  rth 
stiffener. 


For  a  stiffener  located  at  x  = 

■b 


[63] 


Where  Ks  indicates  integration  over  the  cross  sectional  area  of  the  S*-*1 
stiffener.  Substituting  the  expressions  for  <5^. 

We  obtain 

vt  •{  f  {&-*>[*(&)-&)*&'*'<'&) ]  1641 

*  [i  <%)  -  *  d77-])  ^ 

-/V  SASlJArJ, c 

*■  JA,  2 


Now 


fAf  i  ^ 


=  the  moment  of  inertia  of  the  stiff¬ 
ener  about  its  own  neutral  axis  r  _  _ . 

[65] 

where  is  the  distance  from  the 

neutral  plahe  of  the  stiffener  to  the 
centroid  of  the  stiffener  (  =  0) 


Thus 


/  fj(,4  -  £Jg\ 


[66] 


The  terms  containing  Ar  are  the  contribution  from  stretching  and  the  one 
containing  2*r  is  the  contribution  from  bending. 
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Similarly 

v,  j { I  6-*)  i  -5t.  ftp)  1671 
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These  values  of  the  work  done  and  an  appropriate  series  for  the  assumed  de¬ 
flection  can  then  be  used  with  the  variational  principle  or  energy  approach 
to  determine  the  plastic  behavior  ot  the  plate.  A  particularly  appropriate 
series  for  the  rectangular  plate  would  be 


v. 


'y^rrjL. 


H9' 


This  particular  series  can  then  be  substituted  into  the  internal  work  ex¬ 
pression,  V, above  and  the  variational  principle  as  given  by  equation  [14] 
can  be  employed  to  determine  the  set  ot  (in  this  case  mn  takes  the 

place  of  n  in  the  general  equation  given  by  equation  [23]).  In  the  cases 
involving  pulses  whose  duration  are  greater  than  the  fundamental  period  of 
the  plate  we  can  represent  /t*s~  as  a  single  term  series 


✓w-r  *  w;  fa)  ^  2L2r 

L 


3.  An  alternate  approximation  for  large  plastic  deformations 

As  is  seen  in  the  previous  sections  the  elastic  region  is  characterized 
by  bending  only  (eq.  [46]  assumes  bending  stresses  only)  and  the  plastic 
region  is  described  by  equations  involving  both  -bending  and  membrane  de¬ 
formation.  The  large  plastic  deformation  ot  plates  can  be  predicted  by 
using  only  membrane  deformation. The  equation  describing  the  behav¬ 
ior  of  ag  isotropic,  cross  stiffened  or  sandwich  plate  can  be  written  as 
follows: 

.  T*  3 $+  r?  *fC- 

where  ?  [70] 

"Tit.  ~  O*  "Tip.  5  u*-i 

is  the  ultimate  stress  ot  the  material  and  ate  the  equi¬ 

valent  thicknesses  for  stretching  in  the  x  and  y  directions  respectively. 
The  equivalent  thicknesses  for  representative  cases  are  given  in  the  fig¬ 
ure  bej.ow.  _ _ H *3C.  , 


* - - 

Case  A.  Stiffened  Plate 

JL  =  A.  +  A^+  /),,  +  -t~r 


Cas§  B.  Sandwich  opiate 


A.  -  /q»  ±  +  T'W,  A  J 

X  s  ^  -  **/  + 

(A  refers  to  cross  sectional  areas  of  stiffeners) 

Fig.  11  Equivalent  Thicknesses  ot  Plating 
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For  an  isotropic  plate  the  equivalent  thickness  is  the  actual  thickness 
of  the  plate.  The  beha/ior  of  isotropic  membranes  under  dynamic  loading 
was  studied  by  Baker  and  Hoffman30  some  years  ago  and  by  this  author2y 
more  recently.  The  solution  of  equation  [70]  for  a  rectangular  plate 
which  has  its  edges  supported  can  be  written 


_ U) 

-tu  'v'-  ^ 


^JL2r 

L 


[71] 


Substituting  into  the  equation  of  motions  multiplying  both  sides  of  the 
equation  by  j£2l  and  integrating  over  the  area 

of  the  plate  we  find  that  ft )  satisfies  the  following  equation 


h>7*. 


[72] 


The  solution  can  be  written  in  the  same  form  as  the  elastic  solution  given 
in  the  previous  section  of  this  report,  i.e. 

J  *  >  I  ~  =  I  .  -1 


[73] 


^ - ^ 

f*-r) 

where  the  load  distribution  Z3*'*, ir )  was  represented  as  follows: 

prx,. ,*)  =  Gujf)Prt)  [’4] 


and 


- 


Tic  ^  V 
^  a v 


>• 


[75] 


4.  Fitting  the  plate  theories  together 

The  previous  sections  have  given  various  theories  for  describing  the  re¬ 
gimes  of  plate  deformation.  These  theories  can  be  used  in  conjunction 
with  each  other  to  follow  the  complete  behavior  of  the  plate  from  the  elas¬ 
tic  through  the  plastic  region.  The  procedure  is  as  follows: 

a.  Compute  the  elastic  behavior  up  to  the  point  where  yield  starts,  then 

b.  Using  the  initial  conditions  determined  from  the  elastic  regime  and  the 

plastic  relations  given  in  2  or  3,  determine  the  plastic  deformation  un¬ 

til  the  velocity  of  the  plate  reaches  zero  for  the  first  time.  It  will 
be  assumed  that  the  plate  is  "frozen"  at  this  time  in  the  same  manner 
that  Allen5  used  in  the  beam  case.  The  final  plastic  deformation  is 
then  given  by  this  frozen  configuration. 
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C.  Cylindrical  shell  type  structures 

The  theory  of  plastic  deformation  of  cylindrical  shells  under  short  dur¬ 
ation  impulsive  loads  has  been  given  in  a  series  of  reports  by  the  prc*- 
sent  author^' The  general  behavior  was  described  in  the  earlier  sections 
of  this  report.  The  energy  absorbed  or  work  done  by  the  internal  forces 
during  deformation  is  as  follows: 


/  L 
yj  > 


>  /■  i 


7V».)  7+t 


[76] 

ziP)^u 

rfi  \/J  1 


where 


1  VW/T-tf  *V 
—  £~  A  Cj’TT *  Lt 


-VYpfj  *  (&)($)  \  »feyrt)  ls  fir**) 

t  Y  *)•(*)  YZ)Y£)  Y  ££>  W*)  Y£)  Ytjy& )  ^ ) 

+c?)YmifY + 


) 

\ 


In  the  above  relations  a  linear  hardening  material  was  assumed  as  was  done 
in  the  case  of  plates  in  previous  sections  of  this  report  and  it  was  assum¬ 
ed  that  the  deflection  took  on  a  pattern  a«Sc  -J-fx-'jcf)  where  Oc.*  ^/l. 
In  these  equations  a  is  the  radius  of  the  shell,  L  the  length,  t  the  thick¬ 
ness,  0  Poisson's  ratio  . 


For  collapse  type  behavior  where  the  shell  forms  a  hinge  line  straight  aero 
the  center  of  the  shell  on  the  side  facing  the  blast,  the  deformation  pat¬ 
tern  is  described  as  follows: 


O  X  JC  '  <■  4 


**'V*V)  =  a  C4rt~cf  -  \/c*  2-  ^  -  XjL 

For  buckling  the  pattern  is  as  follows 

7Tk.  ‘  C  o<  ef  x.rr 

—  ^  C</a_  f i- TT- cf) 


[78] 


[79] 


Tr  ^  ^  2.77- 


and  finally  for  a  very  short  terms  impulse  where  the  shell  is  loaded  over 
half  the  circumference  and  a  "ring  type"  pattern  results,  we  have 


Fig.  12  Deformation  Pattern  of  Shell  with  Sprayed 
Explosive  Over  Half  the  Circumference 

-k  f'  -  ^  Z^“ *  *  ** 

- 


4  *  "Va, 

*4  <  x.'  x(*/u  +■  efi*7?- L 


[80] 


by  symmetry  we  double  the  energy  obtained  for  the  first  quadrant  (  <fi  <  */'*.  ) 
also  the  energy  for  (Jt  +  */i)<  x.'  <  L.  is  the  same  as  that  for 

X* <  */(_  (also  by  symmetry). 

The  above  functions  have  been  programmed  for  the  computer.  Some  numerical 
results  are  given  in  previous  reports^' and  some  further  results  and  com¬ 
parison  with  experiments  will  be  given  in  the  next  section  of  this  report. 


VII. 


A  more  exact  theory  for  axially  symmetric  elastic-plastic  deformation 
of  cylindrical  shells 


A.  Uniform  isotropic  shells 

The  equation  of  elastic  deformation  of  isotropic  cylindrical  shells 
under  axisymmetric  loading  is 


h  =  thickness,  E  =  elastic  modulus, 


where  O  = 

density,  Pfat  )  is  the  load  distribution. 


[81] 
=  mass 


This  equation  can  be  solved  subject  to  boundary  conditions  at  the 
edges  x  =  0,  JL  and  appropriate  initial  conditons.  For  a  uniform 
shell  made  of  perfectly  plastic  material,  the  shell  becomes  plastic 
when  the  yield  condition  is  reached  at  a  given  locations,  i.e. 
when 


<r»  x  -  ai  o  +  <r<f 1  =  gz  *•  cr*,  j  [82) 

If  this  condition  occurs  at  the  edge  of  a  fixed  beam  then  the  boun¬ 
dary  condition  at  the  edge  is 

w  =  0  (edge  deflection  zero) 

M  =  MQ  (bending  moment  equal  to  plastic  moment) 

If  this  plastic  condition  occurs  at  any  other  section  then  we  must 
resort  to  continuity  conditions  for  beams,  i.o. 


=  n, t 

A 

nt  r  nc 


(deflection  on  the  left  of  the  hinge  =  deflection 
on  the  right  of  the  hinge) 

(left  bending  moment  =  right  bending  moment) 

(left  shear  =  right  shear) 

(bending  moment  equal  to  plastic  moment) 


The  finite  difference  computer  program  developed  several  years  ago 
for  beams'^  can  be  applied  directly  to  this  problem  by  making  changes 
in  the  yield  condition.  It  is  anticipated  that  the  basic  behavior  of 
shells  going  into  "ring  type"  deformation  and  center  hinging  can  be 
ascertained  directly  from  this  theory. 

B.  Uniform  anisotropic  shells  (including  layered  shells) 

The  basic  behavior  of  certain  anisotropic  shells  under  general  non- 
axisymmetric  loading  was  given  by  Flugge^  and  applied  to  some  dynam¬ 
ic  problems  by  this  author.^  The  equation  for  axially  symmetric 
lateral  deformation  can  be  written  down  immediately  from  these  re¬ 
sults.  This  equation  is  the  anisotropic  counterpa-*  .f  equation  [81] 
and  includes  such  factors  as  initial  external  or  !r  rnal  static 

I 


pressure.  It  is  as  follows: 

a+  Ot*r  r  k 


+£s * <*&*&).  £t  f  (izkdkl  1 

«t>, c  /.  pya^  J 

)] 


[85] 


/  - 


&) 


/-  P'/P> 


+  f-f  ^  i±  +  2d.  y.  /p'A  A,  \/-  _  Oji  \  ~1 

L  *v*  *cl  />**<  ^  + i£ jCJEzJe! 

/  -  P‘/0»  J 


x  £>*. 


p* is  the  initial  pressure  in  load/unit  area,  p'  is  the  edge  loading  in 
load/unit  length,^  is  the  mass  per  unit  area,  f*,  it )  is  the  rad¬ 

ial  pressure.  The  appropriate  values  of  the  constants  for  sandwich, 
cross  stiffened,  and  orthotropic  shells  are  given  in  the  previously  cited 
reference.  32  One  case  of  considerable  interest  is  the  layered  shell  which 
is  not  given  specifically  in  the  above  reference-3  but  the  constants  can  be 
written  down  directly  by  using  the  work  of  Flugge31  and  this  reference.32 
First  assume  that  each  layer  of  the  n  layers  is  orthotropic  and  has  a 
stress-strain  law 


0JCw  —  Gcf 

-  Ej*.  A  Ef^  &q) 


[86] 


Then 


where  is  the  elastic  modulus  in  the  x  di¬ 

rection  of  the  nth  layer  and  is  the  thickness  of 
the  layer.  [87] 


1 
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^  •  2^  where  is  the  elastic  modulus  in  the  <p  direction  of 


the  nth  layer 


P.  = 


^  t . 


/r*  =  £  ^  -t*2  ±  Z. 


where  AVs.  is  the  distance  from  the  midsurface  of  the  entire 
shell  to  the  midsurface  of  the  nth  layer 


"  /2. 


Z-  Egf^.  ^"v.  ^ 


-S^  -.5  =  =0  t92l 

and  the  other  constants  are 

JT-  ^  lT=/y  Z3*  '  sj)  J  *  *  1  J  S~~'  [93] 

f  -  I j  *  =  /  >  a-'  jjr**  j  ?*-i  j  ^ ^ 

VTH.Some  simple  applications  and  comparisons  with  experiment 

A.  Beams 

Consider  the  simple  case  of  the  uinform  cantilever  under  blast  from 
pentolite  and  let  us  apply  the  impulse  theory,  neglecting  elastic  ener¬ 
gy.  As  given  earlier  in  the  report  the  energy  flux  (energy  per  unit 
area)  can  be  written  in  terms  of  the  impulse  per  unit  area  as  follows: 

~  r  r  r  s  /  [94] 

Z/*  yti  s  * sj  fie r  g,r  bc*»~ t 


The  value  which  Goodman  obtained  for  the  side  on  impulse  per  unit 
area  was  given  in  a  previous  section  of  this  report  and  its  value  is 


7  _  3  J* f  3  r*illncc  [95] 

fiS*  X  tk  "* 

where  ..  2-  / 

X  -  VT/i2J 

-  atmospheric  pressure  in  ATM 
2  =  **  R  =  distance  from  explosion  to  target 

W  =  weight  of  explosive 

Applying  the  previously  developed  ideas  of  plastic  deformation  of  hing¬ 
ing  cantilever  beams,  the  energy  absorbed  at  the  hinge  is 

V  =.  [96] 
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Assume  a  constant  hinge  moment  M0 


Thus 


V  =  nce> 


where  w  =  width  of  beam 
07  =  yield  stress 
t  =  thickness  of  beam 
A  =  loaded  area  of  beam 

Using  the  appropriate  values  for  the  energy  flux,  the  impulse,  etc.  we 
obtain  _ 

X./J2  J  _  3s*  2  *  ./3  2.J 

d  -  ± J _  [98] 

f*  -  mass  density  of  beam  material 


Consider  an  aluminum  beam  12"  long,  .091"  thick  and  1"  wide.  The  results 
for  the  tip  deflection,  d,  given  by  the  above  relation  and  the  corres¬ 
ponding  experimental  results^"  are  given  in  the  following  table 


Table  I.  Theoretical  and  Experimental  Results  for  Cantilevers 


Charge  Weight 

Distance 

Theoretical  Tip 
Deflection 

Experimental  Tip 
Deflection 

2.  s' 

9" 

h.  m 

3' 

7./" 

7.  3  " 

4' 

+  7" 

2-  7  " 

S' 

3.3" 

f 

3r 

!+  ' 

s  2" 

.  os* 

There  are  sufficient  calculations  shown  in  the 

table  to  illustrate 

trend.  For  small  charges  very  close  to  the  beam,  i.e.  for  very  short 
duration  pulses  where  diffraction  loading  is  of  primary  consideration, 
the  impulse  theory  coupled  with  the  rigid  plastic  assumption  predicts 
good  results.  However  for  large  charges  further  away  where  gust  load¬ 
ing  has  a  substantial  contribution  the  results  predicted  by  impulse  theo¬ 
ry  and  rigid  plastic  behavior  are  not  good.  Therefore  recourse  must  be 
made  to  the  more  complete  theory  such  as  the  finite  difference  approach 
mentioned  earlier. 
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B.  Cylindrical  shells 

1.  Shells  subjected  to  blast  at  various  ranges 

a.  Prediction  using  energy  formula  directly 

We  assume  perfectly  plastic  behavior  of  the  material.  Thus  the  curves 
of  Fig.  8  of  the  earlier  reference^  can  be  used  directly.  For  com¬ 
pleteness  these  curves  are  reproduced  here  in  Fig.  14.  Some  examples 
are  given  below  which  illustrate  the  application  of  the  theory. 

(1)  A  charge  weight  of  389#  of  pentolite  is  exploded  29'  away  from 
a  steel  shell;  the  shell  is  3"0D,  6"  long,  .019"  thick  and  goes  into 
a  collapse  pattern.  The  problem  is  to  predict  the  final  deformation. 
The  experiment  showed  dcfo  1.0 

The  energy  flux  predicted  by  the  formula  given  in  an  earlier  section 
of  this  report  is 


=  /<t.2x/o 


(  A  »/* 


tp  j  r* 


-  3/J-O 


jr  /*. 


The  total  energy  from  the  explosion  will  be 

V  -  Ef.  >  Ap  (Af,  r  fryct+tol  a***  c *  s£*//-L.o) 


\/  -  3 / ST>  X  2  -  JTS,  700  ^ 

The  theory  says  that  this  energy  is  absorbed  in  plastically  deform¬ 
ing  the  shell.  Since  we  are  assuming  a  perfectly  plastic  material, 
assume  that  (JJ  is  the  ultimate  stress,  which  was  found  experimen¬ 
tally  to  be  about  50,000  psi. 


CJ  // 


Thus  we  obtain 

—  \/  SZ  7o  O  *  i/3 

W  S.  _ _ _  -  - 1- - CJ  // 

V  +  L  to3  X  ‘O)**  Xt. 4 

Going  into  the  curve,  a  V  =  11  corresponds  to  a  &°/q  >  >  /  for 
collapse  of  a  shell  with  **/ o  ~  2.  •  This  means  that  the  theory 

predicts  deflections  which  are  much  too  high. 

(2)  Next  consider  a  charge  weight  of  8.4#  at  8'  away.  The  shell  is 
steel,  3"0D,  11.62"  long,  . 019* thick  (collapse  pattern).  The  experi¬ 
ment  showed  d*  /d  ~  <  7 


The  energy  flux  predicted  by  the  formula  is 

Thus  e+  -  r-jxr)  «.  wo  — 

V  -  3/jSvo  * t  ^ 

V  -  (3I.TOC  */.  ,6  1.)  ~ 

Even  in  this  case  a  V  of  3.3  predicts  collapse  deformation  which 
is  much  too  great. 
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(3)  Next  consider  1.1#  at  3'  away  from  the  steel  shell  3"  OD,  11.62"  long, 
.019"  thick  with  a  collapse  pattern.  The  experiment  showed  s,  J 

Ef  -  /,o£  *  to  3  2^2  .  V  -  i.9 


Again  the  prediction  is  too  large. 

Prediction  using  side  on  impulse  to  compute  energy 

Instead  of  working  with  the  energy  formula  supposing  we  use  the  side  on 
impulse-*-'  to  compute  the  energy.  The  energy  will  then  be  given  by 


=  JL 
2. 


cx) a- 


*  2- 

(1)  Consider  the  first  shell  solved  in  the  previous  section.  For  this 
case  the  side  on  impulse,  X  ,  is  about  90  psi  millisec,  thus 


.  S~x  O  *  )  3  _ 


-  1*10 


rt  /»o 
/*»  *• 


V  - 


ok  /•  1—  ^  / 

>3  X.  <5/4  K 


.  IZ+  *■ 

Thus  _ 

\f  ~  ^  *-OQ 

Sit  x /O' 

This  value  of  V  predicts  a  deformation  of  /Q  /  ;  a 

deformation  of  **>/p  ~  / _ was  found  experimentally. 

(2)  Consider  the  second  shell  solved  in  the  previous  section.  For  this 
case  the  side  on  impulse  is  about  25  psi  millisec 

Er  —  '  S  X  ^  -  22,  4-  *  — 

,  y  .  O/Cf  ,h  v- 

\y  —  22.  -f  K  J  A  //.  G  2-  =  7  9-0  ** 

V  ~  J  ^  ^s~b  X.C/C/X.  /.  Sx//'  Cl.  J  ~  >C)  S' 


Coming  over  on  the  curve  of  L/p  -  we  obtain  a 

The  experiment  gave  a  value  of  e^° ~  »  7 _ 


do /&  rx.  ,  £ 


(3)  F 


For  the  third  shell  l  -  14.5  psi  millisec, 

Thus  **+/£>  fer  the  theory  is  .4.  Experiment  shoved 


V 'Z.az  a- 

‘V°/fc  ^  *  5 

(4)  Consider  another  case  of  a  steel  shell  exposed  to  8.4#  explosive  at 
5.8’.  The  shell  is  6"  diameter,  17.5"  long  and  .035"  thick  and  goes 
into  a  collapse  pattern.  The  side  on  impulse  is  about  30  psi  millisec. 

The  V  for  thjLs  case  is  about  .03  which  corresponds  to  a  &*e/o  ,  -4- 

The  experimental  value  of  do/p  was  about  .5. 


(5)  Now  consider  buckling  failure  of  an  aluminum  shell.  The  shell  is  3"0D, 
2"  long,  .008"  thick, subject  to  an  explosion  of  1#  at  a  distance  of  8’. 
We  use  an  experimentally  determined  0^  of  about  15000  psi.  Tha  assumed 
value  of  k  in  the  buckling  formula  [79]  is  k  =  \  .  ‘ 

These  values  were  assumed  for  all  the  subsequent  buckling  calculations 
in  this  reoort. 
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The  side  on  impulse  for  this  case  is  about  8psi  millisec.  The  calculated 
value  of  V  ir  about  .4.  Going  into  the  curve  we  find  that  for  buckling 
the  value  of  is  about  .05.  The  experimental  value  of  *w*/<*.  for 

this  test  was  found  to  be  .08. 

(6)  Next  consider  buckling  of  an  aluminum  shell  3"OD,  5"  long  and  .008"  thick 
subject  to  an  explosion  of  1.1#  at  10'.  The  side  on  impulse  is  about  6 
psi  millisec.  The  calculated  value  of  7  is  about  .3  which  gives  a 
of  about  .07.  The  experimental  value  was  about  .1. 

;(70  Lastly  consider  collapse  of  an  aluminum  shell  3"OD,  9"  long,  and  .042" 
thick  subjected  to  an  explosion  of  115#  at  17'.  The  03  for  this  case 
was  50,000  psi.  The  side  on  impulse  is  about  64  psi  millisec.  The  cal¬ 
culated  v  was  about  .31  giving  a  f  .  The  experimental  value  of 

d*/p  was  about  .6. 


2.  Sprayed  explosive 


In  section  IIB  3  of  this  report  we  discussed  the  behavior  of  shells  under 
very  short  duration  loading  such  as  would  be  obtained  from  sprayed  explo¬ 
sive.  It  has  been  found  experimentally  that  the  fixed  ended  shell  deforrtjs 
into  a  pattern  which  correspv  Is  to  an  abrupt  change  at  the  ends  and  an 
almost  uniform  deformation  along  the  length.  Consider  as  a  first  approxi¬ 
mation  that  end  effects  are  negligible.  The  energy  of  deformation  can 
then  be  written  as  follows:  (assuming  perfectly  plastic  behavior) 


*  *  *  *  iZ*  -  (?)  V.  '¥,>  W  ] 

x2-cp 

Now  for  a  uniform  impulse  over  half  the  cylinder 


Evaluating 


-T  =  T.J>*  =  \/v  ^Vrr^L 

the  integrals  contained  in  \/  above,  we  obtain 


[99] 


[100] 


[101] 


An  aluminum  cylinder  was  tested,  the  parameters  of  which  were 

CL  -  /.  JT  “  =  2>Oj  ooc,  f  si 

s*  -  -/v. -t 

The  experimental  deformation  was  about  *5"  under  a  measured  impulse  of  about 
27  psi  millisec. 

Using  the  above  relation  we  find  that  for  -.2  3  (which  corresponds 

to  V  deflection  in  a  cylinder  with  a  =  1.5") 

1 ~  2  **  >t  ft  sec- 
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Thus  the  theory  predicts  reasonably  well  the  impulse  for  a  given  deflec¬ 
tion  for  this  particular  test. 

C.  Cpnclusions  based  upon  experimental  comparisons 

It  is  shown  that  the  side  on  impulse  criterion  gives  excellent  results 
for  beam  deformation  in  the  case  of  very  small  charges  and  short  distance. 
For  larger  charges  and  distances  we  will  have  to  resort  to  the  more  exact 
elastic-plastic  theory  using  the  differential  equation. 

For  cylindrical  shells  it  seems  that  the  side  on  impulse  gives  excellent 
results  for  rigid-plastic  collapse  and  buckling  behavior  of  all  the  shells 
under  long  or  short  duration  blasts.  When  more  cases  of  sprayed  explosive 
results  become  available,  they  will  be  compared  with  theory. 
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